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INTRODUCTION 


The  cha rac  te  r  i s  t  i  cs  of  a  communication  system 
depend  strongly  on  the  kind  of  the  channel,  e.g.  HP, 
tropospheric  scatter,  microwave,  etc.  Many  of  the  channels 
used  for  spread-spectrum  communications  have  time-varying 
characteristics.  The  channel  considered  in  this  report  is 
the  HP  ionospheric  channel,  and  the  spread  spectrum 
technique  considered  is  frequency  hopping  spread-spectrum. 
The  character ist ics  of  such  a  spread  spectrum  system  are 
determined  and  limited  by  the  multipath  propagation  and 
fading  in  the  HF  channel  as  well  as  by  additive  noise. 

A  signal  transmitted  at  HF  propagates  to  the 
receiver  via  several  paths  determined  by  single  and  multiple 
reflection  conditions  of  the  ionospheric  layers.  Since  the 
propagation  time  associated  with  each  of  the  paths  is 
different  due  to  their  different  optical  path  lengths,  the 
received  signal  suffers  time  spread  distortion  that  degrades 
the  performance  of  the  communication  system.  The  heights  of 
the  reflecting  ionospheric  layers  are  time  variable,  and 
hence  the  signal  is  Doppler  shifted  on  reflection,  thus 
causing  different  frequency  shifts  on  the  various  multipath 
components.  The  presence  of  ionization  irregularities  causes 


variable  diffractive  fadings  on  each  of  the  multipath 
components  and  resulting  fading  of  the  composite  received 
signal.  The  frequency  shift  and  the  fadings  combine  to  cause 
frequency-spread  distortion  that  degrades  the  performance  of 
the  system. 

To  obtain  immunity  to  intentional  interferers  one 
can  design  a  system  which  can  hop  over  a  very  large  band  of 
frequencies.  As  a  result  of  the  large  frequency  separation 
between  hopping  frequencies  we  assume  that  each  hopping 
frequency  propagates  via  statistically  independent  HF 
ionospheric  channels.  For  a  frequency-hopping  communication 
system  operating  via  the  HF  ionospheric  channel  described 
above  we  propose  an  optimum  receiver  and  acquisition  scheme 
whose  characteristics  are  analyzed. 


In  this  work  we  deal  with  the  analysis  of  a 
frequency-hopped  spread-spectum  communication  system  over  a 
non  specular  Rayligh  fading  multipath  channel.  We  assume 
that  the  frequency  separation  between  each  hopping  frequency 
is  large  enough  to  ensure  that  they  propagate  via 
statistically  independent  channels.  This  means  that  each 
transmitted  signal,  with  different  hopping  (carrier) 
frequency,  will  suffer  statistically  independent  time  and 
frequency  spreading. 

Suppose  now  that  we  transmit  a  signal  s^(t)  whose 
hopping  frequency  is  fj.  This  signal  will  propagate  to  the 
receiver  via  several  paths.  The  receiver  receives  the  signal 
Sj(t)  after  the  time  delay  of  seconds.  This  delay  is  the 
propagation  delay  and  depends  on  the  optical  path  lengths 
associated  with  each  path. 

The  signal  s^(t),  with  hopping  frequency  f.., 
propagates  to  the  receiver  via  several  paths  whose  lengths 
are  different  than  the  lengths  in  the  previous  case.  Thus 
the  propagation  delays  t/  and  tj  associated  with  the  signals 
Sj(t)  and  Sj(t)  with  carrier  frequencies  fj  and  f .. 
respectively,  are  different. 

Let  us  now  transmit  s^ft)  with  carrier  frequency 
fj  at  tj  and  tj*  Because  the  heights  of  the  reflecting 
layers  are  time  variable,  the  path  lengths  associated  with 


Sj(t)  at  tj  and  are  different. 

In  a  frequency  hopping  spread  spectrum  system  the 
transmitter  carrier  frequency  is  determined  by  the  output 
sequence  of  a  pseudo  random  generator.  The  number  of 
different  hopping  frequencies  is  determined  by  the  length  of 
the  pseudo  random  sequence.  If  the  pseudo  random  generator 
is  a  maximum  length  linear  feedback  shift  register  (LFSR) 
generator  the  length  is  given  as 

L=2n-1  (1) 

where  n  is  the  number  of  stages  in  the  LFSR  generator.  So, 
each  hopping  frequency  will  be  repeated  after 

T=t2-t1=(2n-l)Th  (2) 

where  T,  is  the  hopping  period.  Therefore  we  can  assume  that 
h 

the  time  delay  between  the  same  hopping  frequencies  is 
sufficiently  large  to  ensure  that  there  is  no  correlation 
between  path  lengths  at  tj  and  t?. 

According  to  these  assumptions  each  hoppinq 
frequency  f.  (l<i<TL)  will  propagate  to  the  receiver  with 
random  delay  time  .  Each  ^  (Ui^L)  is  a  random  variable 
with  probability  density  function 


Pr(t)  * 

V. 


2  At 


Tnd-At<t<Tnd+At 


(3) 


elsewhere 


where  T  .  is  the  nominal  or  mean  delay  time,  and  At  is  the 
nd 


4 


maximum  deviation  from  T 


Since  any  two  frequencies 


nd  • 

propagate  via  statistically  independent  channels ,  the  random 
variables  and  are  uncorrelated  i.e. 

E{(  Z[  -If)  (  tj  i  *  j  (4) 

Also  the  random  variables  at  t^  and  ate  uncorrelated 
i.e. 

Ef  r  *|(tj)-  ^(tj)i  r  T[ft2)-  fyt2n }  =o 

The  nominal  delay  time  Tn^ ,  in  the  general  case, 

is  a  function  of  frequency.  Pig.  1  gives  typical  curves  for 

Tnd  versus  frequency  for  three  different  paths. 

In  Fig.  2,a,b,c,d  and  e  are  sketches  of 

measurement  results  of  the  time  delay  of  a  signal 

propagating  via  ionospheric  multipath  channel.  As  we  can  see 

from  Fig.  2,c,d  and  e  it  is  possible  to  identify  typical 

curves  for  nominal  delay  time  Tnd.  But  it  is  impossible  to 

recognize  the  typical  curves  of  the  nominal  delay  time  from 

Fig.  2, a  and  b.  The  reason  for  this  is  that  the  path  lengths 

are  random  variables  and  therefore  the  delay  time  associated 

with  each  path  is  a  random  variable  whose  mean  value  is  Tnf? , 

which  of  course  is  a  function  of  frequency. 

Suppose  now  that  the  frequency  hopping  rate  is  of 

the  same  order  as  At  (maximum  deviation  of  the  delay  time 

from  T  . ) .  In  such  a  case  there  will  be  significant 
n  d 

overlapping  between  two,  three  or  more  hopping  frequencies 
at  the  receiver  side.  With  such  overlapping  of  multiple 


adjacent  hopping  frequencies,  fading,  and  white  additive 
Gaussian  noise,  standard  FH  receivers  will  neither  acquire 
or  track  the  received  signal,  nor  will  such  a  receiver 
detect  the  data  with  a  reasonable  error  rate.  The  problem  of 
optimum  detection  and  acquisition  is  solved  in  this  report. 


DELAY  TIME 


t 


h 

t 


REVIEW  OF  PRIOR  WORK 


The  synchronization  problem  is  extremely  important 
in  any  spread-spectrum  communication  system.  The  initial 
synchronization  or  acquisition  is  the  most  difficult  of  all. 
Specific  synchronization  requirements  depend  mostly  on  the 
intended  application,  and  often,  the  synchronizing  scheme 
has  to  be  designed  for  the  worst  case  condition.  Suitable 
measures  of  system  performance  are  the  mean  time  of 
acquisition  and  the  probability  of  successfully  acquiring  an 
anticipated  spread-spectrum  signal.  The  scheme  to  be 
employed  generally  depends  on  the  application,  the  amount  of 
time  allowed,  and  the  amount  of  uncertainty  involved.  The 
probability  of  acquisition  is  a  suitable  performance  measure 
in  situations  where  the  transmitter  continuously  emits  the 
spread  spectrum  signal  that  the  receiver  must  acquire.  The 
mean  time  to  acquire  is  a  good  criterion  for  burst 
spread-spectrum  systems. 

Many  techniques  for  achieving  synchronization  have 
evolved.  Below,  we  review  some  of  the  frequency-hopping 
acquisition  schemes. 


-  11  - 
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Matched  filter  scheme 


A  near-optimum  Matched  Filter  for  detection  of  the 
frequency-hopped  signal  sequence  is  given  in  Fig.  3a.  Each 
arm  of  the  filter  contains  a  mixer,  a  bandpass  filter,  a 
square  law  detector  and  delay  line.  Assume  that  the 
spread-spectrum  signal  hops  over  M  distinct  frequencies,  and 
that  the  frequency  hopping  sequence  is  f1,f2,...,f^  which 
repeats  itself.  The  acquisition  scheme  then  consists  of  M 
mixers  each  followed  by  a  square  law  detector  and  delay 
line.  The  delay  lines  are  inserted  so  that  when  the  correct 
sequence  appears,  the  maximum  of  the  voltages  Vj,V2,...,VM 
will  occur  at  the  same  time.  Therefore,  the  output  voltage 
of  the  adder  will  exceed  the  threshold  level  with  high 
probability,  indicating  synchronization  of  the  receiver  to 
the  signal. 


Stepped  serial  search  scheme 


The  stepped  serial  search  acquisition  scheme  is 
given  in  Fig.  3b.  This  scheme  is  used  for  acquisition  of 
long  pseudo  random  codes.  In  a  stepped  serial  acquisition 
scheme  the  locally  generated  pseudo  random  frequency  hopped 
signal  is  correlated  with  the  incoming  frequency-hopped 
signal.  When  the  local  hopping  is  aligned  with  that  of  the 
incoming  signal,  the  input  to  the  envelope  detector  is 


ideally  a  sinusoid  at  the  intermediate  frequency,  and  the 
output  of  the  integrator  will  exceed  the  threshold  with  high 
probability. 

Both  schemes  are  analyzed  in  f 1 1 r  and  their 
performances  are  given.  The  performance  measure  is  the 
probability  that,  at  the  end  of  the  examination  interval, 
the  desired  signal  code  sequence  is  not  detected.  The  miss 
probabilities  versus  signal-to-noise  ratio  for  a  Ricean 
fading  channel  are  given. 

These  two  acquisition  schemes  are  basic 
acquisition  schemes  for  frequency-hopped  spread  spectrum 
signal.  A  two  level  threshold  scheme  is  also  analyzed  in  [11# 
and  [2].  In  [3]  and  [41  the  sequential  estimation  technique 
is  described. 

The  commonality  amtheong  these  analyses  of 
acquisition  schemes  is  that  they  use  the  same  model  for  the 
HF  channel.  It  is  assumed  that  the  propagation  delay  time  is 
constant  for  all  hopping  frequencies  and  that  the  spreading 
time  is  not  a  function  of  time  and  frequency.  This 
assumption  is  acceptable  for  slow  frequency-hopped  signal 
where  the  overlapping  {due  to  difference  in  delay  tine) 
between  two  adjacent  frequencies  can  be  neglected. 


SUMMARY  OP  RESEARCH  COMPLETED 


Tn  this  work  we  propose  an  acquisition  scheme  for 
a  frequency-hopped  SS  signal  received  over  a  Rayligh  fading 
multipath  channel,  whose  delay  time  versus  frequency  is 
given  in  Fig. 2.  We  consider  relatively  fast  FH  with  a 
hopping  rate  of  1  KHz,  and  127  different  hopping 
frequencies.  We  assume  that  the  probability  density  function 
of  the  channel  delay  time  is  given  by  (3),  with  the  maximum 
deviation  from  Tnd  of  1  msec.  The  received  signal  is 
embedded  in  white  Gaussian  noise.  While  the  above  problem 
appears  to  be  specific  the  approach  given  below  is  general 
and  can  be  generalized  further. 

In  this  case,  ignoring  the  spreading  effect,  there 
will  be  a  possible  overlaping  between  only  three  adjacent 
frequencies.  For  such  a  case  we  propose  an  optimum  receiver 
and  acquisition  scheme  whose  characteristics  are  analyzed 
mathematicaly  and  by  computer  simulation. 

We  characterize  the  performance  of  the  acquisition 


scheme  by  the  probability  of  detection  for  given  probability 
of  false  acquisition  when  the  input  signal  is  white  Gaussian 
noise  only.  The  probability  of  detection  is  the  probability 
that,  at  the  end  of  the  observation  period,  the  desired 


signal  code  sequence  is  detected. 
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RESEARCH 


The  system  we  consider  is  shown  in  Fig. 4.  The 
transmitter  consists  of  a  frequency  synthesizer  whose  output 
frequency  is  determined  by  n+m  bits.  The  m-bits  are  data 
bits  and  the  n-bits  are  bits  generated  by  the  LFSR 
generator.  Actually  this  is  a  case  of  an  m-ary  FSK 
frequency-hopped  SS  system. 

When  there  is  no  input  data  the  output  frequency 
is  determined  by  the  LFSR  generator  only.  This  case  is 
assumed  when  acquisition  is  to  be  performed. 

Receiver 


Let  us  for  a  moment  assume  that  the  channel  shown 
in  Fig. 4.  consists  of  a  variable  delay  line  only,  and  that 
the  nominal  delay  time  Tnd  is  not  a  function  of  frequency. 
Actually  we  assume  that  the  curve  Tnd  versus  frequency  is 
known  (Fig.  1). 

The  transmitted  signal  is  given  by 

Sj(t)=A  cos{  (w^Wj )  t]  l<i<L  ;  (6) 

For  this  case  the  received  signal  for  At  =  0  is  shown  in 
Fig. 5(a)  and  the  received  signal  for  -Th<At<Th  is  shown  in 
Fig. 5(b).  In  this  figure  A  means  that  a  frequency,  say  fj 
can  arrive  advanced  by  a  time  as  much  as  T,  and  D  means  that 


INCOMING 
FH  SIGNAL 


Fig.  6 


it  can  be  delayed  by  as  much  as  . 

We  can  see  from  Fig. 5(b)  that  there  is  possible 
overlapping  between  three  adjacent  frequencies  only,  in  the 
time  interval  3T^  to  4Th,  for  example,  there  is  possible 
overlapping  between  f._2,  fi-i  anrf  f i  *  So  at  the  same  time 
we  can  receive  three  signals  with  three  different  carrier 
frequencies.  In  order  to  receive  complete  signal  without 
this  overlapping  we  need  three  receivers  Rj,  R?  and  R3  whose 
input  circuits  are  adjusted  to  f  .  _  3  ,  fj  and  f.+3 
recpect ively .  Assuming  perfect  synchronization  i.e.  assuming 
that  we  know  exact  time  of  arrival,  for  each  frequency,  for 
At* 0  (no  variable  delay),  the  receiver  which  will  detect  all 
signals  without  overlapping  is  shown  in  Fig.#>. 

Synthesizers  1,2  and  3  will  hop  to  the  new 
frequency  each  3Th  seconds  following  the  sequence  shown  in 
Fig. 5(c) . 

Acquisition  scheme 

The  basic  function  of  an  acquisition  circuit  is  to 
bring  the  incoming  PN  sequence  and  the  locally  generated 
sequence  into  sufficiently  close  alignment  so  that  the 
difference  is  within  the  pull-in  range  of  the 
fine-acquisition  (tracking)  loop. 


Let  us  again  assume  that  the  channel  shown  in 
Fig. 4  can  be  modeled  by  a  variable  delay  line.  We  also 
assume  that  during  the  initial  acquisition  no  data  is 


transmitted,  i.e.  the  transmitted  signal  is 


s . ( t) =A  cosfw.t) 


0<t<!Th,  Ui<L 


(7) 


In  this  case,  the  receiver  does  not  know  the  exact  time  of 
arrival  of  the  signal,  but  the  receiver  knows  that  the 
signal  has  to  arrive  in  the  time  interval  from  n  to 


(2  -1)T^+2T^  seconds.  Also  the  receiver  knows  the  sequence 


of  the  transmitted  frequencies.  Thus  if  we  are  able  to 
determine  the  time  of  arrival  *or  any  frequency  for  at=0, 
the  synchronization  will  be  perfect. 

So,  we  have  a  case  of  detection  of  a  signal  with 
unknown  time  of  arrival,  amplitide  and  phase.  The  optimum 
receiver  for  this  case  is  shown  in  Fig. 7a.  There  are  three 
quadrature  receivers  because  there  is  possible  overlapping 
between  three  adjacent  frequencies  only.  For  the  transmitted 
sequence  given  in  Fig. 7b,  the  outputs  and  are 
shown  in  Fig. 7c.  The  delay  lines  of  T^  and  2T^  will  ensure 


that,  for  a  noiseless  channel,  the  maxima  of  and 


Ui+1  will  occur  in  the  interval  3Th  to  BT^ . 


Acquisition  strategy 


During  the  time  interval  T= (2n-l ) T^+2T^  the  output 
signals  U j _ j • U ^  and  ui+i  are  observed  continuously  and 
compared  with  the  threshold  TR.  If  any  output  signal  axceeds 
the  threshold  TR  the  maximum  is  detected  and  the  time  of  its 


occurance  T.  (i=l,2,3)  is  determined. 


If  only  one  maximum  is  detected  we  assume  that 
this  maximum  is  due  to  the  noise  only,  and  we  assume  that 
acquisition  has  not  occurred. 

In  the  case  when  two  maxima  are  detected  and  when 
the  time  distance  between  them  is  greater  then  3T^  we  again 
declare  that  acquisition  has  not  occurred.  If  this  distance 
is  less  then  3T.  we  calculate  the  acquisition  time  as 


ACQUISITION  TIME  = 


Ti  +  Ti 


When  all  three  maxima  are  detected  and  when  the 

time  distance  between  T.  and  T.  (i  =  j)  is  less  then 

1  ,max  3  ,m ax 

3T.  we  calculate  the  acquisition  time  as 


ACQUISITION  TIME  » 


Tx  +  T?  +  T3 


If  Ti  mav  “  T-i  is  greater  then  3T,  and  if  there  is  no 
i  finax  3  f m l n  n 

time  distance  between  any  two  maxima  less  then  3T^ ,  we 
declare  that  acquisition  has  not  occurred.  If  there  is  such 
a  distance  the  acquisition  time  is  calculated  as  in  the  case 
with  two  maxima. 


False-alarm  probabilit 


The  probability  of  false  acquisition  Pfacn  1S  the 
probability  that  acquisition  will  occur  if  the  input  signal 
to  the  receiver  of  Fiq.7a  is  the  noise  only.  Before  we 


we  have  to  find  the  false  alarm  probability 


calculate  Pfacn 

for  a  system  that  observed  the  output  q(t)  of  a  quadrature 
receiver  over  an  interval  of  time  O^t^T  which  is  much  longer 
than  the  duration  of  the  signal.  A  false  alarm  occurs  in 
this  system  if  the  output  q(t)  exceeds  the  threshold  level 
TR  at  any  time  during  the  observation  interval,  when  the 
input  to  the  system  consists  of  noise  alone.  Let  us  denote 
the  probability  of  this  event  by  P^a.  Stated  otherwise  this 
probability  is  given  by 

Pfg  =  1  -  P (T)  (10) 

where 


P (T)  =  Pr { q ( t) <TR ,  0{t«T}  (11) 

is  the  probability  that  the  receiver  output  q(t)  does  not 
appear  above  the  threshold  TR  at  all  during  the  interval 
O^t^T.  We  can  assume  that  at  t=0,  q(0)=0,  i.e.,  that  the 
receiver  output  will  start  from  zero.  In  this  case  for  TR>0, 
P(0)=1  and  P (oo)  =0 . 

The  negative  derivative 

IP 

p(t)  =  -  -  0<t<  oo  (17) 

it 

is  the  so  called  "first  passage-time  probability  density 
function";  p(t)dt  is  the  probability  that  the  output  q(t) 
crosses  the  threshold  TR  from  below  for  the  first  time  in 
the  interval  t  to  t+dt.  The  problem  is  to  calculate  the 


density  function  p(t).  This  problem  can  be  solved  for  only  a 
few  types  of  stochastic  processes. 


The  first  passage  time  probability  density 
function  can  be  found  for  the  case  of  a  one  dimensional 
Markov  process. 

Jt  can  be  shown  that  the  output  of  the  quadrature 
receiver  given  in  Fig. 8  is  a  Markov  process  whose 
transitional  probability  density  function  is 


p(qt/q) 


ft 

&*(!-/<*) 


•  expf- 


2  6za-/t2)  °  n-fi2’)£z 


(13) 


2 

where  o  is  the  variance  of  the  noise  at  the  output  of  each 
RC-filter,  and  /£=exp(-f  S'1)  •  It  can  be  shown  (71,(81  that  for 
normalized  variance  C  =1,  (13)  satisfies  a  Fokker-Planck 
differential  equation 


dj> 

bt 


\  1  A 

-2—  {K1(q)p(q)  }  + - (K2  (q)  p ( q)  > 

3  q  2  q 


with 


(14) 


] 

Kj  (q)  - - q  K?(q)  =  2  (15) 

q 

The  Fokker-Pianck  equation  can  be  used  to  find  the 
probability  density  function  p(TR,t/qQ)  of  the  first  time  t 
that  the  envelope  q(t)  crosses  the  threshold  level  q=TR 
given  the  value  q  at  t=0.  The  first  passage  problem  can  be 


solved  by  the  method  of  Siegert  rsi  .  According  to  Siegert, 
the  Laplace  transform  of  the  f i rst-passage-t ime  p.d.f. 
p(TR,t/q)  is 


2*P<TR,t/q  )  } 


W( -  ;  1 ;  ) 

2  2 

S  TR2 

M( ;lf  ) 

2  2 


j 


V 


where  M(a;b;z)  is  a  confluent  hypergeometric  function  and 
according  to  [61  is  given  by 


az  a(a+l)  z.  a(a+l)  (a+n-1)  z 

M  (a ;  b;  z)  * -  +  -  — —  +•  •  •+  -  — II-  + 

b  b  (b+1 )  2!  b(b+l)  (b+n-1)  n! 


In  order  to  find  the  probability  density  function  itself  we 
must  find  the  inverse  Laplace  transform  of  (16).  By  the 
usual  technique  one  can  obtain  a  series  of  the  form 


p(TR,t/qQ)  *^~Anexp(5pt) 


where  the  A  are  the  residues  of  (16)  at  the  roots  S  of 
n  n 


;)  -  0 


For  values  of  (TR  /2)>1  the  largest  zero  of 
(19)  is  greater  than  -1  and  such  that 


^ 


where  S.,  is  the  second  largest  zero.  Tn  this  case  the 


27 


equation  (18)  can  be  approximated  by 

p(TR,t/qQ>  *  AjexpfSjt)  (21) 

From  (18)  and  (12)  we  can  find  the  probability 
that  the  output  of  a  quadrature  receiver  from  Fig.P  will 
not  cross  the  threshold  level  TR  during  the  interval  0<t<T, 
will  be 


P(T) 


T) 


(22) 


n  n 


or  using  the  approximation  (21) 


P(T) 


exp(SjT) 


(23) 


and  finally  using  (10)  ,  the  probability  of  false  alarm  is  given  as 


Pfa  -  1  +  — —  exp(SjT) 
S! 


(2*) 


Usually  we  need  to  know  the  value  of  the  threshold 
for  a  given  probability  of  false  alarm. 


Calculation  of  the  threshold  TR  for  a  given  P^3 


Using  (22)  we  can  express  the  probability  of  false 


alarm  as 


- 1 


exp ( ST) 


(2M 


where  An  are  the  residues  of  (IS)  at  the  roots  Sn  of  (19). 
In  this  case  the  probability  of  false  alarm  is  given  and  the 
period  of  observation  T  is  known.  We  also  assume  that  qQ=0. 
It  follows  that  in  order  to  find  TR  we  have  to  determine 
M  ( S/2 ; 1 ; TR j/2  )  so  that  the  roots  of  this  confluent 
hypergeometric  function  and  the  residues  of  (IS)  at  these 
roots  satisfy  (25)  for  a  given  P^g  and  T. 

This  problem  has  been  solved  numerically  since  it 
is  impossible  to  find  a  solution  in  closed  form. 

For  different  values  of  Pfg  and  for  T  =  129  msec 
the  normalized  threshold 


TR  = 
n 


has  been  calculated  and  is  given  in  Table  1.  and  is  plotted 
in  Fig. 9. 


Table  1 


Pf  a 

TR 

n 

0.9 

3.55 

0.5 

4.02 

0.25 

4.25 

0.1 

4.52 

0.075 

4.50 

0.05 

4.70 

0.025 

4.85 

0.01 

5.07 

0.0075 

5.1?  . 

0.005 

5.21 

0.0025 

5.35 

0.001 

5.54 
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Probability  of  false  acquisition  P 


racn 


False  acquisition  will  occur  if  the  time  delay 
between  any  two  maxima  exceeding  the  threshold  is  less  then 
or  equal  to  3T^,  when  the  input  signal  is  noise  alone.  Let 
us'  denote  by  Wj(t^)  (i*l,2,3)  the  maximum  due  to  noise  only 
occurring  at  and  with  M..(t2)  f j=l , 2 , 3 ; i»j)  the  maximum 
due  to  the  noise  only  occurring  at  t2*  The  probability  of 
false  acquisition  can  then  be  calculated  as 

I 

1 

I  Pfacn“  3P^  |ti“t2|  ^3Th}  PfW>TR}  P{M2(t2)>TR)  (27) 

|  where 

P{M j  f t j ) >TR)  =  PfMj ft2)>TR}  =  Pfg  (2P) 

is  actually  the  probability  of  false  alarm.  For  TR=0  the 
probability  of  false  alarm  is  Pfa=l  and  (27)  becomes 


D  — 

facn 


3P<  IV'J  °v 


(29) 


Since  the  probability  of  false  acquisition  P, 

1 3cn 

is  usually  given,  we  have  to  set  up  the  threshold  TR  so  that 

this  given  condition  is  satisfied.  In  order  to  calculate  the 

threshold  TR  for  given  P-_„  we  first  calculate  P,  from 

racn  fa 

(27) .  The  result  is 


1  3P{  |tl“t2| <3Th> 

When  we  have  P^a  for  a  given  Pfacn  »the  normalized  threshold 
can  be  found  from  Fig. 9,  and  by  using  (25)  the  true  value  of 
the  threshold  can  be  calculated. 

Probability  of  detection 

The  input  signal  rft)  to  the  quadrature  receiver 
of  Fig. 8,  is  given  by 

r  ft)  *  sft  -  ?)  +  nft)  f?l) 


where  nft)  is  Gaussian  white  noise  and 


s/t  -?*)  *  A  cosfwft  -T)  +B) 


0<t<Th 


where  ©  is  a  random  phase  with  uniform  distribution  between 
0  and  27 T  and  T"  i  s  a  random  time  of  arrival  with 
probability  density  function 


p(T) 


0<r«!T 


elsewhere 


According  to  the  theory  Of  optimum  detection  of  a 
signal  with  random  time  of  arrival  T91  we  have  to  observe 
the  output  of  a  quadrature  receiver  (Fig.R)  from  0  to  T,  to 


detect  the  maximum  and  the  time  of  its  occurrence.  The 
maximum  is  to  be  compared  with  the  threshold  TR ,  and  if  that 
value  exceeds  the  threshold  we  say  that  signal  is  present 
and  the  time  at  which  the  maximum  has  occurred  is  the  time 
of  arrival.  The  value  of  the  threshold  is  to  be  chosen 
according  to  the  given  value  of  false  alarm  probability  Pfa- 
In  order  to  make  the  calculation  possible  we 
assume  that  the  maximum  due  to  the  noise  and  signal  occurs 
at  f  +Th  i.e.  seconds  after  the  signal  appears  at  the 
input  of  the  quadrature  receiver.  The  value  of  this  maximum 
is  a  random  variable  with  a  Ricean  probability  density 
function 


p(q) 


exp(- 


2  .2 
q  +  A 


■)i  > 


where  Aq  is  the  maximum  value  of  the  signal  at  the  output 
when  there  is  no  noise  at  the  input;  IT  is  the  variance  of 
the  noise  at  the  output  of  the  RC  filter. 

The  probability  that  the  signal  will  be  detected 
is  given  as 

Pd  ®  P{sig+noise>TR  and  noise<sig+noise}  (35) 

The  probability  that  the  niose  will  not  reach  the  value  of 
the  envelope  at  t*  ?'+Tjl  during  the  time  interval  n<t4T  is 


exp(S_T) 


q  -  is  the  value  of  the  signal  plus  noise  at  t»  +T^  and  is 
a  random  variable  with  Rice  density  function  p(q) .  Averaging 
over  all  values  for  q  greater  than  the  threshold  TR  we  get 
the  probability  of  detection  as 


or 


P  (q) p(q) dq 


T* 


(37) 


oo 


p 


d 


exp(SnT)p(q)dq 


(38) 


r« 

In  this  case  A  and  S_  are  functions  of  q.  Equation  (38)  has 
n  n 

been  solved  numerically.  For  different  values  of  Pfa  the 
probability  of  detection  P^  has  been  calculated  and  the 
calculated  results  have  been  compared  with  the  results 
obtained  by  computer  simulation.  In  the  simulation  we  have 
assumed  that  the  signal  and  its  time  of  arrival  are  detected 
if  the  maximum  due  to  both  the  signal  and  noise  is  in  the 
interval 


T  +7^/2  <  t  <  r+  3Th/2  (39) 

is  greater  than  any  maximum  due  to  the  noise  only.  Simulated 
and  calculated  results  are  given  in  Fig. 10, 11,12, and  13 


Probability  of  acquisition 


The  probability  of  acquisition  PgCq  can  be  shown 


to  be 


P _ -  3(1-P„)P*  +  p! 


(40) 


Calculated  and  simulated  results  are  given  in 


Pig.  1*5, 17  and  18. 

Probability  of  false  acquisition,  PfJe 


False  acquisition  will  occur  if  at  least  two 
detected  maxima  are  due  to  the  noise  only  when  the  signal  is 
present,  and  their  time  distance  is  less  than  3Th .  In  order 
to  calculate  Pfac  we  have  to  find  the  probability  that  the 
noise  is  greater  than  threshold  and  greater  than  signal  plus 
noise  at  t*?'+Th.  We  can  write  this  as  follows 


P 


fas 


P{sig+noise>TR  and  noise*sig+noise J 
P{sig+noise<TR  and  noise>TR> 


+ 


(41) 


This  probability  can  be  calculated  as 

OP  TR 


Pf as*  ^  fl-Pfq)  )p(q)<Jq  +  Pfa  ^  p (q) dq 


(42) 


TR  0 

and  finally  after  few  manipulations  we  have 

TR 


Pfas*(Pfa"1)  J  P(q>dc*  ♦  1  ~  ^  P(q)p(q)dq 


(43) 


0  TR 

Now  we  calculate  the  probability  of  false 


acquisition  Pfac  as 


Pfac  "  (Pfas>  °‘139 


(44) 


where  0.139  is  Pfacn  for  TR*0. 

Calculated  and  simulated  results  for  Pfas  and  Pfac 
are  given  in  Fig.  14,15  and  Fig.  16,17,18  respectively. 


Fig. 10  Probability  of  detection  versus  signal  to 
noise  ratio  for  Pfa“  1 
-  calculated  results 


simulated  results 


PROBABILITY  OF  DETECTION  IN  PERCENT 


PROBABILITY  OF  DETECTION  IN  PERCENT 


vu 


Fig. 12  Probability  of  detection  versus  signal  to 


noise  ratio  for  P^a«0.268 


i 


-  calculated  results 

-  simulated  results 


<dB> 


SIGNAL-TO-NOISE  RATIO  E./N 

n  o 

Probability  of  acquisition  Pacq  and  probability 
of  false  acquisition  Pjac  Pfacn“  0.01 

-  calculated  results 

-  simulated  results 


IN  PERCENT 


SIGNAL-TO-NOISE  RATIO  E.  /N  (dB) 

n  o 

Probability  of  acquisition  Pacg  and  probability 
of  false  acquisition  Pfac  for  given  probability 
of  false  alarm  P^a*l. 

-  calculated  results 


simulated  results 


CONCLUSION 


The  maximum  probability  of  false  acquisition  Pfacn 

is  for  TR»0,  and  this  value  is  0.139,  and  is  same  for  all 

values  of  signal  to  noise  ratio. 

As  we  can  see  from  the  diagrams  showing  for 

given  probability  of  false  alarm  Pja»  the  difference  between 

simulated  and  calculated  results  is  acceptable  in  the  entire 

range  of  signal  to  noise  ratio.  The  difference  is  larger  for 

lower  signal  to  noise  ratio.  The  reason  for  this  is  that  for 

small  signal  to  noise  ratio  there  is  more  probability  that 

the  maximum  due  to  the  signal4-noise  will  be  shifted  to  the 

left  or  to  the  right  from  t*2-'+Th.  For  high  signal  to  noise 

ratio  the  difference  is  almost  negligible. 

For  the  same  reason,  the  difference  between 

calculated  and  simulated  results  for  p  _  is  larger  for 

acg  3 

lower  signal  to  noise  ratio. 
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FUTURE  WORE 


In  this  work,  we  have  analyzed  an  acquisition 
scheme  for  a  frequency-hopped  SS  communication  system  over 
an  HF  channel  whose  delay  time  is  given  by  (3),  neglecting 
the  effect  of  multipath.  To  get  the  complete  characteristics 
of  the  system,  the  multipath  spreading  has  to  be  included  in 


analyses.  For  this  case,  the  calculations  of  P.,  P  and 

d  acq 


Pjac  are  impossible  and  therefore  a  computer  simulation  has 
to  be  performed.  The  consequence  of  multipath  is  to  increase 
the  acquisition  time. 

Also  the  tracking  problem  for  such  a  channel  has 
to  be  solved.  Since  in  this  system  overlapping  between  two, 
three  or  more  adjacent  frequencies  is  possible,  the  tracking 
problem  is  much  more  difficult  than  in  the  case  of  constant 
delay  time. 


With  the  scheme  we  have  proposed  the  minimum 
acquisition  time  is  approximately  equal  to  the  period  of 
pseudorandom  sequence.  For  very  long  sequences  the 
acquisition  time  can  be  unacceptable.  Therefore,  future  work 
should  include  the  development  of  rapid  acquisition  schemes. 
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